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Abstract 

A time-marching finite difference code, 
XTRAN3S, that solves the three-dimensional 
transonic small perturbation equation for flow 
over isolated wings has recently been developed. 
During initial applications of the program, 
problems were encountered in the prediction of 
unsteady forces. The use of a revised grid and 
force calculation scheme improved those predic- 
tions. Comparisons are made between predicted 
and experimental pressure data for a rectangular 
supercritical wing. Comparisons of steady and 
unsteady data at Mm = D.7nn show good agree- 
ment between calculated and experimental values. 
A comparison of steady data at Mm = n.825 
shows poor agreement between calculations and 
experiment. Program difficulties have been 
encountered with swept and tapered configura- 
tions. 


Nomenclature 
R ' wing span, inches 

p 

c nondimensional local chord,— 

c n section normal force coefficient 

c r reference chord, inches 

C local chord, inches 

C m wing pitching moment coefficient 

C n wing normal force coefficient 

C p pressure coefficient, 

AC p lifting pressure coefficient (difference 

between lower and upper surface pressure 

coefficients) 

f oscillation frequency, Hz 

F airfoil contour 

c to 

k reduced frequency, 

k nondimensional time scaling 

Moo freestream Mach number 

p local static pressure 

poo freestream static pressure 

q» freestream dynamic pressure* 

kVT 

t nondimensional time, — — 4 

c 

T time, seconds r 

Vco freestream velocity, inches)second 

x,y,z cartesian coordinates nondjmensionalized 

h V c r 

X,Y,Z cartesian coordinates, inthgs 
a angle of attack, degrees > 

Aa pitch oscillation amplitude*, degrees 


Y ratio of specific heats 

£,n,C nondimensional transformed coordinates, 

(eq. 8) 

rf fractional span, 

a) frequency of oscillation, radians/second 

<J> phase angle between pressure and wing 

motion, degrees (positive for pressure 
leading motion) 

$ nondimensional perturbation velocity 

potential 

velocity potential function, 

3 = VmC r (x+<!>) 

[ ] denotes jump in quantity across a 

d-i s conti nuity 

Subscri pts 

LE leading edge 

TF trailing edge 

o initial or mean value 


Introduction 

In recent years there has been significant 
development of finite difference methods for 
solving the unsteady transonic flow equations. 
These methods have direct application to flutter 
and other aeroelastic problems which tend to be 
most severe in the transonic range where linear- 
ized flow analyses are inaccurate. For two- 
dimensional flows, methods based on the tran- 
sonic small perturbation (TSP) equation are well 
developed, fine of the most widely used programs 
is LTRAN21-2 which is used to solve the low 
frequency time-dependent TSP equation. The 
LTRAN2 code has been extended in its frequency 
ranged- 4 and has been further extended to the 
complete time-dependent equation by including 
the 4>tt term. 5 These methods have been used 
extensively for aeroelastic analyses (for 
example, see refs. 5-7 and their references). 

The time-dependent solution algorithm has 
been extended to treat the three-dimensional TSP 
equation. The resulting code, XTRAN3S, 
treats isolated wings and includes unsteady 
aerodynamic and aeroelastic analysis capability. 
Results are given for the low frequency equation 
in refs, in and 11 and for the full frequency 
equation in refs. R, 9 and 12-15. The addition 
of boundary layer effects has also been demon- 
strated.! 4 

XTRAN3S has been implemented on the CPC 
CYRER. 203 at the NASA Langley Research Center 
(LaRC). This paper describes initial investiga- 
tions that have been conducted to develop and 
validate the CYBER 203 version of the program. 
Results of the program are evaluated by compari- 
son with other methods and experiment. 
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First, a brief review of the TSP equation 
and boundary conditions solved in XTRAN3S is 
presented. Next, assessments of and changes in 
the finite difference grid and the method of 
calculating aerodynamic forces are discussed. 
Calculations of steady and unsteady pressures 
are shown for a rectangular wing with a super- 
critical airfoil section that was recently 
teste ?R ™ t ^ ie T rans °nic Dynamics Tunnel at 
LaRC.15 Comparisons of XTRAN3S predictions 

with the corresponding test data for a limited 
number of cases are presented. Program 

instabilities encountered for swept and tapered 
configurations are discussed. 


XTRAN3S Program Description 

The XTRAN3S program can he used to solve 
static and dynamic problems for rigid and 

elastic wings in transonic flow. A brief 

summary of the governing equation and boundary 
conditions given in refs. 9 and in is presented 
here. 

A time accurate alternating direction 

implicit (ADI) finite difference scheme is used 
to solve the three-dimensional unsteady 

transonic small perturbation (TSP) equation 

-(M^V 2 »£&x>t + tD-'&’V A *X + R? y) x 

+ (* y + <V y ) y + (* z ) z = n ( 1 , 


A = B = C = 0 


(5) 

The boundary conditions imposed 

upon the 

flow field are 


Far upstream: 

4 = n 

(6a) 

Far downstream: 

4 x +k$t = n 

(6b) 

Far above and below: 

4 Z = n 

(6c) 

Far spanwise and wing root: 4y = n 

(6d) 

Trailing vortex wake: 

G 

II 

r ~N 

(fie) 

[*X 

+ k* t ] = 0 

(fif) 


The airfoil flow tangency condition is 


$ z = F x + kF t ; z = n+ * X LE £ x £ X TE (7) 

The far field boundary conditions (6c) and (fid ) 
are equivalent to treating the computational 
boundaries as solid walls. 

The wing is mapped into a rectangular car- 
tesian coordinate system using the shearing 
transformation 

x-x 

* («a) 

n = y (8b) 


The coefficients A, R, and C depend upon 
the assumptions used in the development of the 
TSP equation. If the unsteady TSP equation is 
derived to be analogous to the steady, three- 
dimensional equation developed by Lomax et 

al. 16 , the coefficients are defined as 

A = - j (y+ 1)M F (2a) 

R = \ (2b) 

c = -(Y-l)M^ (2c) 


These are referred to in ref. in as the "NASA/ 
AMES" coefficients. Deriving the TSP equation 
from the mass conservation equation as proposed 
by van der Vooren et al. 17 leads to the “NLR" 
coefficients 


A 

R 


- 1 CM?-y}m f )m^ 




(3a) 

(3b) 


S = z (Rc ) 

This maps any swept or tapered wing into a 
rectangular planform in the computational 

domain. For swept and/or tapered planforms this 
leads to a computational grid that is skewed in 
the physical domain. 

The program is run at LaRC on a GDC CYRER 
203 using a computational grid with 60x20x40 
points in the £,n, and c directions. The CYRER 
203 is capable of scalar or vector arithmetic. 
The scalar version of the XTRAN3S program 
requires 1.9 central processing unit (CPU) 

seconds per time step. Vectorizing the Ani 
S-sweep and a matrix manipulation routine 

reduces the required time to 1.2 CPU seconds per 
time step. Further savings in CPU time are 

possible by reorganizing the Ani algorithm. For 
a rigid rectangular wing with a 6% parabolic arc 
airfoil section, the solution converges to a 
steady state answer in 720 time steps and con- 
verges to a harmonic solution in B cycles (1440 
time steps). 


C = 



(3c) 


The "classical" coefficients based upon an 
analogy to the two-dimensional steady transonic 
von Karman-Guderl ey equation are 

A=-|(Y+1)M l • (4a) 

R = C = fl (4b) 

Finally, the coefficients for tfu4 ^linear equa- 
tion are 


The XTRAN3S program can be used to solve 
for aeroelastic motions of the wing through a 
coupling of the aerodynamic solution and the 
structural equations of motion. A linear model 
of the structure is used and a numerical 
integration of the structural equations is 
performed for each time step in the aerodynamic 
calculation. The structural deformations are 
represented using generalized coordinates 
defined in terms of structural mode shapes. 

XTRAN3S Program Verification 
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During initial applications of XTRAN3S, 
calculated forces exhibited spurious oscilla- 



tions during convergence to a steady state solu- 
tion. The cause of the oscillations was traced 
to the computational grid. I 8 In ref. 1 8, an 
analysis of the accuracy of linear unsteady 
forces calculated on several two-dimensional 
grids was conducted. The study was carried out 
using a two-dimensional TSP finite difference 
program, XTRAN2L. Included in the grids analy- 
zed was the default x-z grid from XTRAN3S. A 
revised x-z grid was described in ref. 18 which 
improved the accuracy of the unsteady calcula- 
tions. 

The first step of the program verification 
was to determine if accurate solutions are cal- 
culated when the code is run in a two- 
dimensional mode. This was done hy repeating 
the grid study of ref. 18, using XTRAN3S with 
the default and revised x-z grids. Comparisons 
are made between calculated and theoretical 
values of the linear unsteady forces to deter- 
mine program accuracy. The grid study was then 
repeated for a three-dimensional case to verify 
that the results of the two-dimensional analysis 
apply in three-dimensions. 


Unsteady Force Calculation 

In order to determine the accuracy of 
predicted forces, the linear unsteady aero- 
dynamic forces on a flat plate were computed for 
a wide range of reduced frequencies. The pre- 
dictions are compared to accurate solutions of 
the subsonic kernel function lifting-surface 
theory. 19-2H Unsteady forces are typically 
determined by calculating several cycles of 
forced harmonic oscillations and performing a 
Fourier analysis on the last cycle of oscilla- 
tion. The complex fundamental component of the 
Fourier analysis provides the force estimate. 
A second approach is the indicial response 
method in which a Fourier transform is applied 
to the force response resulting from a step 
change in a given mode of motion. A transient 
pulse technique, 7 * 18 a variation of the indi- 
cial response method, was used to determine the 
unsteady forces for the present study. 

In the transient pulse technique, a 
smoothly varying exponentially shaped pulse for 
a given mode of motion, (e.g. angle of attack) 
is applied to the wing. The unsteady calcula- 
tion is carried out until all transients in the 
forces have decayed. A fast Fourier transform 
(FFT) is then applied to the input and force 
response time histories. The force transform is 
then divided by the transform of the input to 
obtain the unsteady aerodynamic forces for a 
wide range of reduced frequencies. The pulse 
used is described in ref. 18 and results in 
accurate force estimates for reduced frequencies 
of k < 2.n. 


Two-nimensional Analysis 

Grid Effects . - The linear analysis was 
performed for a two-dimensional flat plate at 
Me = n. 85 n, a 0 = n.no, v«/c r = sn.n 

sec- 1 , k = l.n and At = n.122718. The time 


step, At, was chosen to give a frequency incre- 
ment between calculated unsteady forces of 
Ak = n.02B. To use XTRAN3S for linear problems, 
the coefficients defined in eq. ( R) were used. 
Two-dimensional results were obtained by defin- 
ing the spanwise computational grid to consist 
of five points, four on the wing and the fifth 
at the computational boundary. The two- 
dimensional linear results are compared to an 
accurate subsonic kernel function analysis. 19 

The first case analyzed is the calculation 
of unsteady forces using the default x-z grid 
described in refs. 9 and 18. A similar analysis 
was performed using a two-dimensional finite 
difference program and is detailed in ref. 18. 
The comparison between unsteady forces calcu- 
lated using XTRAN3S and those obtained from the 
kernel function analysis is shown in Fig. 1. As 
in ref. 18, forces calculated using the default 
grid exhibit very poor agreement with kernel 
function calculations over a wide range of 
reduced frequencies. The oscillations exhibited 
in the calculated unsteady forces are due to 
reflections of waves that emanate from the 
wing. For k < 0.3 the oscillations are caused 
by reflections from the z boundary, which is 

treated as a solid wall (eq. fic). For k > 0.3 
the oscillations are caused by reflections from 
internal grid points due to the spacing of the 
computational grid in the z direction. 

A revised x-z grid was proposed in ref. 

18. This grid is designed to improve the calcu- 
lation of unsteady forces by reducing the 
reflection of waves. The revised grid includes 
modified spacing in the z direction to minimize 
internal reflections, and the computational 
boundaries were moved further away from the wing 
to reduce the frequency and magnitude of the 

boundary reflections. The unsteady forces 
calculated using the revised x-z grid are shown 
in Fig. 2. The results obtained using the 

revised grid show much better agreement with 
those from the kernel function method for a wide 
range of reduced frequencies. The internal and 
boundary reflections have been significantly 
reduced, resulting in a much better comparison 
of forces for k _< n.R. The forces are 
inaccurate for k > n.S due to the continued 
existence of some internal reflections. 

Force Calculation Techniques . For reduced 
frequencies In the range n _< k _< n.s, an offset 
exists between unsteady forces calculated using 
XTRAN3S and the kernel function method as shown 
in Fig. 2. The cause of this discrepancy was 
traced to the procedure used to calculate the 
forces. Forces are calculated by integrating 
pressures along a chord, which for the normal 
force coefficient is 

c n = /AC p (9) 

In small perturbation theory the pressure 

coefficient is defined as 


The value of $£ is computed at a grid point 
using linear differencing upstream and down- 
stream of the grid point and averaging the two 
values. Integration of the pressures to obtain 



forces is accomplished using Simpson's rule. 
The combination of first order differencing with 
second order integration can lead to an error in 
the force calculation. 

The force calculation can be improved by 
integrating the potential instead of the pres- 
sures as is done in LTRAN2-NLR. 3 For the 
normal force coefficient, combining eqs. (9) and 
(10) gives 

c n = - | - ?V fM t d? (11) 

where A$ is the difference between lower and 
upper surface potentials. The first integral 

can be exactly evaluated, and is equal to 

- 1 - m le> u*) 

Since a$le = n.n, the normal force coeffi- 
cient is 


c n = - \ * a *TE - ?k / A4 t (13) 

Instead of differentiating to calculate A$r 
and integrating to calculate the forces, the 
potential jump at the airfoil trailing edge, 
A$te» is used in obtaining the forces. 

A comparison of kernel function calcula- 
tions and the unsteady forces calculated by 
XTRAN3S with the revised grid and alternate 
force calculation scheme is shown in Fig. 3. 
The revised force calculation has eliminated the 
offset in unsteady forces between XTRAN3S and 
kernel function calculations. For a flat plate, 
the forces show very good agreement for reduced 
frequencies up to k = 0.6. 


Three-Dimensional Analysis 

The linear analysis was performed for a 
flat plate rectangular wing at ^ 0.860, 
a 0 = 0 . 00 , Voo/Cp = 60.0 sec- 1 , k = 1.0 
and At = 0.122718. The rectangular wing has a 
panel aspect ratio of 3.30. As in the two- 
dimensional analysis, the coefficients defined 
in eq. (6) were used. The spanwise computational 
grid is the same as that described in ref. 9, 
consisting of 16 equi-spaced points on the wing 
with an additional point defining the wing tip 
and 4 grid points outboard of the wing. The 
three-dimensional XTRAN3S results are compared 
to an accurate calculation using the three- 
dimensional subsonic kernel function.? 0 

Unsteady forces were calculated using 
the default x-z grid on the rectangular flat 
plate using the revised force calculation 
method. Comparison between the calculated 
unsteady force and kernel function predictions 
is shown in Fig. 4. As in the two-dimensional 
case, forces calculated using the XTRAN3S 
default x-z grid exhibit very poor agreement 
with kernel function calculations. The kernel 
function results for k > 1.6 are erratic due to 
a program limit in the numerical quadrature 
used. 


Using the revised x-z grid with the same y 
grid as above leads to the results shown in 

Fig. 6. The XTRAN3S calculated forces agree 

well with the kernel function forces for 
k <0.6. The cause of the low frequency 

discrepancy is unknown, but the cause of the 
discrepancy at high frequencies is the same as 
in the two-dimensional analysis - reflections of 
disturbances within the computational domain. 
The analysis was performed using several 

spanwise grid spacings, and the calculations 
were found to be only slightly influenced by the 
spanwise spacing. 


Comparisons With Experiment 

In this section, comparisons of XTRAN3S 
calculations and experimental data obtained from 
a rectangular wing recently tested in the 
Langley Transonic Dynamics Tunnel 1R are pre- 
sented. The wing has a panel aspect ratio of 
2.D and a U% thick supercritical airfoil sec- 
tion^ Pressure measurements were made at four 
spanwise stations over a broad range of Mach 
number and angle of attack for both static pitch 
and oscillations in pitch about 0.46 chord. The 
test Reynolds number was 4.0 million based on 
chord and the test medium was Freon (y = 1.131). 
Further description of the tests and the'modei 
are given in ref. 16. All calculations were 
made using measured airfoil ordinates and for 
measured root angle of attack and Mach number. 
The revised x-z grid was used for all subsequent 
calculations. 


Steady Results 

The steady state calculations are compared 
with measured data for a 0 = 2. no and for 
M» - 0.70D and 0.826 in Figs. 6 and 7 respec- 
tively. Calculations were made using the AMES 
(eq. 2) and NLR (eq. 3) coefficients in the TSP 

equation. For VL = D.70D (Fig. 6) there is a 

weak shock near the inboard leading edge. The 

agreement with the data is generally good over 
the mid portion of the chord but with some 

deviation near the nose and the lower trailing 
edge regions. There is little difference in the 
results obtained using the NI.R or the AMES 

coefficients except near the shock. The NLR 

coefficients give a somewhat stronger shock for 
this case. Some of the lack of agreement in the 
nose area may be a result of using a relatively 
coarse grid near the nose. 

For Moo = 0.826 (Fig. 7), the shock is 

further aft on the inboard portion of the wing 
and approaches the leading edge at the tip, 

showing a large three-dimensional effect. The 
correlation of the computed results with the 

experimental data shows trends similar to the 
results for M* = 0.7DD, but with signficantly 
poorer agreement. The NLR coefficients give a 
shock that is significantly stronger and located 
further aft than that predicted using the AMES 
coefficients. 

A sample two-dimensional calculation has- 
been made for this airfoil using a non- 
conservative full potential program 21 with 
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and without a boundary layer. The results for 
Mo© = n.825 and ctp = 0° are shown in Fig. 
8. There is a qualitative correspondence of the 
trend shown with and without a boundary layer 
(Fig. 8) to the correlation between the inviscid 
XTRAN3S calculations and the measured data at 
Moo = 0.825 (Fig. 7). This suggests that vis- 
cous effects are significant for this supercri- 
tical wing. The boundary layer addition to 
XTRAN3S14 has not yet been implemented at 
Langley. It would appear that further effort in 
viscous calculations is necessary. 

It appears that better correlation between 
the inviscid full potential (Fig. 8) and XTRAN3S 
calculations (Fig. 7) is obtained when the NLR 
coefficients are used. At the inboard span 
station, the results from XTRAN3S with the NLR 
coefficients show a shock at the trailing edge, 
which agrees with the two-dimensional full 
potential prediction. 


Unsteady Results 


Hue to the poor agreement in steady state 
pressures betwen experiment and XTRAN3S calcula- 
tions at Moo = 8.825, unsteady data comparisons 
between experiment and calculations are shown 
only for Mo, = 0.7nn. At this Mach number the 
best comparison of steady state pressures 
.between experiment and XTRAN3S calculations is 
obtained using the NLR coefficients (see Fig. 
8). Using the NLR coefficients gives a stronger 
shock on the inboard region of the wing, which 
comes closer to matching the experimental data. 
Rased on this observation, the NLR coefficients 
were used for all unsteady analysis at 
M» = 0.780. 

Unsteady results for M*, = 0.700, 

a o = 2.n° f Act - + 10 , k = n.178 and 

f = ID Hz are shown in Fig. 9 in the form of 
magnitude and phase of the lifting pressure 
coefficient. Corresponding results are shown 
from XTRAN3S using the NLR coefficients and from 
linear lifting surface theory (RHOIV)22. The 
pressure amplitudes calculated by XTRAN3S are in 
good agreement with the experimental data over 
most of the wing with some over prediction in 
the inboard leading edge region and under 
prediction over the outboard portion of the 
wing. The linear theory results are in good 
agreement with experiment except near the 
leading edge where transonic effects are 
evident. As for the pressure amplitude, the 
XTRAN3S results for phase show an over 
prediction in the inboard leading edge region. 
An examination of the experimental data reveals 
that the large change in the phase data occurs 
over the trailing edge portion of the wing's 
lower surface and is probably caused by viscous 
effects. The linear theory results for phase 
are in good agreement with experiment over the 
forward part of the airfoil. The hest agreement 
between XTRAN3S calculations and experimental 
data occurs in the mid-span region where the 
dynamic shock is not over predicted and tip 
effects are not pronounced. 

Similar results are shown in Fig. in for 
k = 0.356 and f = 20 Hz. The comparison with 
XTRAN3S is improved for this higher reduced 
frequency, particularly in the phase near the 


trailing edge. The linear theory pressure 
amplitude prediction is not as good as that for 
f = in Hz, under predicting the pressure near 
the leading edge and over predicting the pres- 
sure near the trailing edge. 


Numerical Stability 

As previously discussed, an An I scheme is 
used in XTRAN3S. For rectangular wings the pro- 
gram exhibits good numerical stability, becoming 
unstable only for large At. Several attempts to 
analyze wings with moderate sweep and taper have 
been unsuccessful due to a numerical instability 
which appears to be independent of At. The 
instability has been encountered for wings with 
taper ratios of approximately 0.3 and leading 
edge sweep from 290-360. Pianforms for 

which the instability has been encountered 
include those of the F-5, AV-8B and nAST ARW-2 
wings. Successful calculations have been made 
for the LANN wingld w hi ch has a taper ratio of 
0.4 and leading edge sweep of 27 . 50 . Thus 

only wings of low to medium sweep and taper have 
been successfully treated. 


, . i j i 1 jf ib m u u Liear 

at this time. It may be related to the highly 
skewed physical computational grid which results 
from the shearing transformation (eq. 8a) or 

to how some of the terms are treated in the 
finite difference scheme. For a wing with 
moderate sweep and taper, a typical x-y plane of 
the grid in physical coordinates is shown in 

Fig. 11. The grid is highly skewed in the up- 

stream and downstream regions, and the grid 
points are closely spaced outboard of the wing. 
This can possibly lead to computational diffi- 

culties. The problem is currently being 

^ssed, but until it is resolved program 
applications are limited to wings of low sweep 
and taper. K 


Concluding Remarks 

A comparison of unsteady forces calculated 
using the XTRAN3S program and accurate theoret- 
ical methods for a flat plate airfoil demon- 
s v n r J ® inaccuracies in the forces predicted by 
XTRAN3S. The cause was found to be the reflec- 
tion of waves from the computational boundaries 
and internal grid points. A new grid signifi- 
cantly reduced errors caused by these reflec- 
tions. A revised scheme for calculating 
unsteady forces improved the agreement of 
calculated and accurate kernel function values. 
A comparison of calculated unsteady forces and 
those predicted by a kernel function method for 
a flat plate rectangular wing demonstrated good 
agreement using the revised grid and force cal- 
culation scheme. 

Comparison of calculations for a rectangu- 
lar supercritical wing with wind tunnel data 
showed mixed agreement. Comparisons of steady 
and unsteady pressures for K» = 0.700 showed 
good agreement while comparisons of steady pres- 
sures for Mo> = 0.825 showed poor agreement. 
Steady state calculations at H, = 0.825 showed 
a large effect of the equation coefficients used. 
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Comparing the experimental data with calcula- 
tions pointed out the need for viscous effects 
In an XTRAN3S-type code. 

Experience with of other configurations 
using the XTRAN3S program demonstrated a numeri- 
cal instability. The instability occurred for 
wings with moderate sweep and taper. Further 
improvements are needed to extend the applica- 
bility of the program to these cases. 


l^Myers, M. R. # Guruswamy, p._ and 
Goorjian, P. M. : Flutter Analysis for a Trans- 

port Wing Using XTRAN35. AIAA Paper 83-nq?? 
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Fig. 1 - Unsteady forces calculated for a flat 
plate airfoil using default XTRAN3S x-z 
grid, M» = .850. 
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Fig. 2 - Unsteady forces calculated for a flat 
plate airfoil using revised XTRAN3S x-z 
grid, H» = .850. 



Fig. 3 - Unsteady forces calculated for a flat 
plate airfoil using revised XTRAN3S x-z 
grid and revised force calculation 
scheme, M», = .850. 



Fig. 4 - Unsteady forces calculated for a flat 
plate rectangular wing using default 
XTRAN3S grid, K» = .850. 
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Fig. 5 - Unsteady forces calculated for a flat 
plate rectangular wing using revised 
XTRAN3S grid, M* = .850. 


7 



*. .1 

U^PER LOWER 

V\ 0 measured 

‘*fr- XTRAN3S (EQ. 2 COEFF.) 

XTRAN3S (EQ. 3 COEFF.) 



Fig. 6 - Comparison of measured and calculated steady state pressure 
distributions for rectangular supercritical wing, M. = 0.700, 
“o = 2°. 
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Fig. 8 - Comparison of invlscid and viscous full 
potential calculations of steady state 
pressure distribution for supercritical 
airfoil, Moo = 0.825, a 0 = 0°. 



x/c 



x/c 


XTRAN3S 

RHOIV 



X/C 


XTRAN3S 

O measured 



Fig. 9 - Comparison of measured and calculated 
unsteady pressures for rectangular 
supercritical wing, = 0.700, 
a 0 = 2°, Act = TlO, k = 0.178 
(Ref. 15, Fig. 14). 
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Fig* in - Comparison of measured and calculated 
unsteady pressures for rectangular 
supercritical wing, = n.7nn, 
a 0 = 2°, Ao = *1°, k = 0.356. 
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Fig. 11 - Typical x-y plane of computational 

grid for wing with moderate sweep and 
taper. 
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